Spitzer Space Telescope evidence in NGC6791: no super-mass-loss 
at super-solar metallicity to explain helium white dwarfs?^ 

Jacco Th. van Loon^, Martha L. Boyer'^, Iain McDonald^ 

66 ■ ABSTRACT 

O 

o. 

^ ■ We use archival Spitzer Space Telescope photometry of the old, super-solar 

p • metallicity massive open cluster NGC 6791 to look for evidence of enhanced mass 

"^ ■ loss, which has been postulated to explain the optical luminosity function and 

Q\ • low white dwarf masses in this benchmark cluster. We find a conspicuous lack 

■ of evidence for prolificacy of circumstellar dust production that would have been 



^ 



c^ 



expected to accompany such mass loss. We also construct the optical and infrared 



Oh' luminosity functions, and demonstrate that these fully agree with theoretical 

O ' expectations. We thus conclude that there is no evidence for the mass loss of 

^ • super-solar metallicity red giants to be sufficiently high that they can avoid the 

^, ' helium flash at the tip of the red giant branch. 

^ . Subject headings: stars: evolution — stars: luminosity function, mass function 

0\ , — stars: mass loss — stars: Population II — open clusters and associations: 

V/-^ , individual (NGC 6791) — infrared: stars 

O. 

95 . 1. Introduction 

o 

. !^ I Low-mass stars in galactic globular clusters lose ~ 20% of their mass during the ascent of 
>< : thefirst red giant branch (RGB; |Dupree et aJ B; l Origha et al. lB; [ McDonald fc van l J 



20071 ). offeri ng a partial explanation for the blue extent of the core-helium burning horizon- 



tal branch (iRood Ill973l ). It has been argued recently that in the old (8 Gyr), super-solar 
metallicity (Z ~ 0.05) massive open cluster NGC 6791, at least 40% of RGB stars experi- 
ence higher rates of mass loss; this would lead to their departure from the RGB before the 
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helium flash, resulting in underma ssive helium white dwarfs fjCastellani fc Castellani Ill993l : 



Hansen I l2005l : iKalirai et al.l 120071 ). If true, this would have important implications for the 



white dwarf populations produced in metal-rich stellar systems, and for the mass return 
into the interstellar medium of such syst e ms — which include giant elliptical galaxies and 



the bulges of spiral galaxies. iBedin et al.l (120081 ) . on the other hand, propose an alternative 
explanation for the white dwarf properties in NGC 6791 without the need for enhanced RGB 
mass loss. 

We investigate the evidence for super-mass-loss on the RGB of NGC 67 91 by looking for 



circu r nstellar dust that m ay form in the winds from cool RGB stars (e.g., iGehrz fc Woolf 
197ll : lOriglia et al.l 120071 ) in archival Spitzer Space Telescope infrared (IR) images, and by 
comparing the IR and optical luminosity functions with model predictions. 



2. Spitzer photometry of NGC 6791 

The observations of NGC 6791 presented her e were intended 



'or James Webb Space 
Telescope photometric and astrometric calibration (iDiaz-Miller 1120071): we obtained the data 



Werner et al. 



2004 ) p ublic archive. The 



from the Spitzer Space Telescope (hereafter Spitzer, 
observations were made with the Spitzer Infrared Array Camera (IRAC; iFazio et al.ll2004l ) 
on 2005 June 14 UT. The image maps cover ^ 20' x 15', with 3.6 and 5.8 fim mosaics offset 
to the northwest from the cluster center by ~ 8' and by ~ 3' at 4.5 and 8 /im. The region 
around the cluster is displayed in Fig. [H High dynamic range mode had been implemented 
to prevent saturation of brighter sources (1 s exposures) while maintaining a high signal-to- 
noise on fainter sources (26.8 s exposures), and observations had been made with an eight 
position cycling dither pattern to build redundancy against outliers and artefacts. 

Raw data were processed with Spitzer Science Center (SSC) pipeline version S14.0.0. We 
corrected the Basic Calibrated Data images for array dist ortions and mosaicked the r esulting 
images with the SSC Legacy MOPEX software vl6.2.1 (JMakovoz &: Marleau I l2005l ) . which 
includes background-matching to minimize pixel offsets in overlapping regions of the mosaics 
and eliminates cosmic rays and other outliers. 

We carried out poi nt-spread function photometry using the DAOphot II photometry 
package (IStetson 1 119871 ). Sources at least 4cr above the background were extracted. A 
sharpness cut-off helped eliminate extended sources and remaining outliers. A pixel-p hase- 
dependent photometric correction was applied to the 3.6 /xm fluxes (JReach et al.ll2005l ). The 
flux estimates are color corrected using a 5000 K blackbody, as specifled in the IRAC Data 
Handbook, version 3.0. Zero magnitudes were taken from the same reference, and flux 
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uncertainties include both the DAOphot errors and cahbration errors fJReach et al.ll2005l ). 
False star tests indicate that the catalogue is complete to at least 90% in the core of the 
cluster down to [3.6] = 17.5 mag and [8] = 16.5 mag. 



3. Cluster membership 



Proper motions have been measured for rt iany stars in and around NGC 6791 (JMonet et al. 



20031 : IZacharias et al.l l2004l : iDias et al.l |2006| ) . We designate a star a proper rn otion non- 



member if the absolut e proper motion i s > 20 mas yr~^ in the lMonet et al.l (120031 ) catalogue, 
> 14 mas vr~ ^ in the IZacharias et al.l (120041 ) catalogue, or its membership probability was 
set zero in the iDias et al.l (120061 ) catalogue. Fig. [2] displays the resulting selection of possible 
members (top row) and proper motio n non-membe r s (^bo ttom row). The figur e illust rates 
the choice of selection criteria for the iMonet et al.l (120031 ) and IZacharias et al.l (120041 ) cat- 
alogues, which are based on the clearly defined clump of stars with small proper motions 
(pos sible members) c ompared to the much more widely distributed foreground population. 
The iDias et al.l (|2006l ) selection criteria include a positional element. Note that if a star was 
excluded as a member by any one of these works it is labelled a non-member in all three 
panels of Fig. [2l 

The proper motion non-members are identified with crosses in the Spitzer [3.6], [3.6]- 
[8] color- magnitude diagram ( Fig. [3]). This eliminates for instance the very brightest star 
(#2240 in lstetson eral]l2003h and the bright very red star at [3.6]=9.7 mag, [3.6]-[8]=2.4 



mag (#2349 in IStetson et al.ll2003l ). Objects in the faint, red corner of the diagram are 
background galaxies; slightly extended objects with similar colors are prominent on the 
Spitzer images. As expected for galaxies, these objects are spread uniformly across the 
images. 

To eliminate field stars with no or ambiguous proper motion information, a further 
criterion is applied on the basis of the optical V, B-V diagram (Fig. HJ using data from 
Stetson et al.ll2003l ): the main sequence, RGB and red clump (around V = 14.6 mag, B—V ~ 
1.35 mag) of NGC 6791 are well separated from the bulk of the foreground main-sequence 
stars. In our subsequent analysis we keep all remaining possible cluster members that have 
proper motion information, photometry in all Spitzer IRAC bands and the V-band, and that 
are located within 0.1° projected distance from the cluster centre (290.23, +37.77). 
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Evidence for circumstellar dust emission 



In the absence of longer-wavelength data, [3.6]-[8] is the most sensitive Spitzer color 
to the presence of circumstellar dust. Dust causes a reddening of this color, firstly, due to 
emission at 8 /xm, and secondly, due to extinction at 3.6 fim if the dust envelope is very 
optically thick. The RGB in NGC6791 is surprisingly narrow in [3.6]-[8] color (Fig. [3]) down 
to at least [3.6] ~ 13.5 mag, i.e. more than four magnitudes below the RGB tip. Compared 
to the distribution mirrored around the peak (Fig. [5]), at [3.6] — [8] = —0.07 mag, there 
are 14 out of 72 stars (19 ± 5%) that are redder than expected from a distribution which is 
symmetric in color. Not only is this but a small fraction, the colors are also only marginally 
red, with A([3.6] - [8]) < 0.1 mag for all but one with A([3.6] - [8]) ~ 0.22 mag. The latter 
is at the faint end of this RGB selection; indeed, there is no sign of the brightest stars be ing 
more likely to show an IR excess as it is observed, e.g., in wCentauri (JBoyer et al.ll2008l ). 



The fraction of dusty stars along the upper 4. 5 mag of the RGB in NGC6791 may be 
compared with that in the globular cluster 47Tuc (lOriglia et al.l 120071 ). There, ~ 340 stars 
are counted in the central part of 47Tuc, amongst which 90 are dusty (26%). In the outskirts 
of 47Tuc, there are ~ 160 stars, amongst which 9 are dusty (6%). These statistics may be 
affected by blending, in the central part, and field star contamination, in the outskirts. 
Nonetheless, there is no evidence for an abnormally high fraction of dusty stars in NGC 6791 
despite an order of magnitude higher metal content than in 47Tuc. 



Following iGroenewegen I (120061 ) . A([3.6] — [8]) = 0.1 mag corresponds to a mass-loss 
rate of 1 — 2 x 10~® Mq yr~^ for an M-type star with aluminium-oxide or silicate dust and a 
luminosity L = 3, 000 Lq, wind speed f = 10 km s~^, and dust-to-gas mass ratio tp = 0.005. 
In the extreme case of an MlO-type star with pure aluminium-oxide dust the corresponding 
mass-loss rate is nearly 5 x 10~^ M0 yr~^, but to our kn owledge no such R GB star exists. The 
mass-loss rate scales approximately as M oc vy/L/ip ( van Loon II2OOOI ). For a dust-driven 
wind one expects the dust-to-gas ratio in the metal-r ich stars of NGC 679 1 to be higher by 
a factor ~ 2.5, and the wind speed by a factor ~ 1.6 (jMarshall et al.ll2004j ). The luminosity 
range along this part of the RGB is ~ 10^~^ Lq, which lowers the Groenewegen value for 
M by up to a factor 5. If the winds are as fast as sometimes seen in chr omospherical lines 
of warmer RGB stars (JMauas et al.l l2006l : iMcDonald fc van Loon 1 120071 ) then M could be 
higher. We thus settle on a conservative estimate of M<10^^ M© yr~^. The mass-loss 
rate may still be higher if the dust-to-gas ratio is lower than assumed — in spite of the 
high abundance of condensab le materia l . How ever, our estimate compares f avorably with 
the empirical relationships of iReimers I (119771 ) and ISchroder fc Cuntz I (120051 ) , which yield 
M ~ 10"^ and a few 10"^ Mq yr'^ respectively, for a typical RGB star of 1 Mq, 4000 K 
and 500 Lq. This suggests that the basic assumptions we made are not unreasonable. 
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If such mass loss is sustained over 19% (the observed fracti on of dusty stars) of the 
time it takes to evolve along this part of the RGB, t ~ 8 x 10^ yr (JMarigo et al.ll2008l ). then 
they will lose AM < 0.2 Mq — insufficient to avoid the helium flash. These stars were born 



with a mass M 



initial 



1.1 Mq, and w ould need to s hed A MinipH eH — 0-7 M^ t o pro duce 



white dwarfs of the masses measured by iKalirai et al.l (120071 ) . Both iKalirai et al.l (120071 ) and 
Hansen I (120051 ) point out that not all RGB stars may avoid the helium flash. Given the 
above estimates of M, fraction of dusty stars and RGB lifetime, the Spitzer data suggest 
that no more than Mt/AMimpiied x 19% or <22% of RGB stars in NGC6791 can sustain 
the required heavy mass loss. 



There is a possibility that we have missed a star at the very tip of the RGB losing mass 
at a high rate, but only briefly before it leaves the RGB. Indeed, dust-accompanied mass 
loss, a t rates as high as 10~ ^ M^ yr~^, is seen to happen predomin antly near the tip of the 
RGB ( van Loon et al.ll2006l : lOriglia et al.l 120071 : iBoyer et al.ll2008l ). It is unclear, however, 
whether such star has sufficient time to lose enough mass to avoid the helium flash. In 
any case, despite the important effect on the remaining mantle mass and its subsequent 
appearance on the horizontal branch, the white dwarf descendant would have a mass that 
is indistinguishable from that of a CO core following th e helium flash, as the luminosi ty on 
the RGB is a direct measure of the helium core's mass (jCastellani fc Castellani Ill993l ). 



5. The RGB luminosity function 

To explain white dwarf masses 0.1 M© lower than the core mass a t the tip of the RGB, 
the lu minosity function must be depleted over 0.6 dex in log L (Table 2 in lCastellani fc Castellani 
19931 ). corresponding to 1.5 mag. We demonstrate below that such depletion in NGC6791 
is not corroborated by the measurements. 

We constructed cumulative luminosity functions (LFs), in both the Spitzer IKAC bands 
and the V-band (Fig. [H]). The IR and V-band LFs were aligned such that the very obvious 
bump due to the red clump coincides — this required a shift towards brighter values of the 
V-band LF by 3 mag. This difference, and the much steeper bright end to the V-band LF 
compared to the IR LFs, are due mainly to the bolometric corrections to the V-band (BCy) 
of the increasingly cooler RGB stars as t hey approach the R GB tip (a small contribution is 
due to the extinction, E{B - r) ~ 0.14; IKalirai et al1l2007l ). 



The observed LFs compare favorably with model LFs from the Padua group (JMarigo et al. 



20081 ). displayed in the lower panel of Fig. [6] The models, displayed for a distance modulus 
of 13.0 mag, show the same shape and location of the red clump and main sequence (the 
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rise in the V-band LF at the faint end of the scale). The models, which do not include mass 
loss along the RGB, reproduce the observed difference at the bright end between the IR and 
V-band LFs rather well. The most metal-rich Padua model is with Z = 0.03 not quite as 
metal-rich as NGC6791 {Z = 0.05), and if compared to the Z = 0.02 Padua model it is 
evident that a Z = 0.05 model LF would be truncated at even fainter magnitudes, in close 
agreement with the observed LF. This is corroborate d by the good fi t to the optical colour- 



magnitude dia g ram of a Z = 0.06 is ochrone using the iReimers I (119771 ) mass-loss prescription 



dClaret 1120071 ). iKalirai et all (120071 ) find that the V-band LF of NGC6791 is depleted with 



respect to the (sub-)solar metallicity clusters Berkeley 17 and M67, over the t op magnitude 



of th e RGB. They do not, however, consider that BCy also differs by > 1 mag (JMarigo et al. 



20081 ). nor the effects from contamination by field stars. 



With regard to the IR LF, which more accurately measures the bolometric output of 
these stars, the difference with the models can be entirely explained by stochastic effects. 
The NGC6791 cluster is not populous enough to harbor many stars near the RGB tip; while 
there are few tip-RGB stars present now, it is conceivable that one or two more may have 
been present several million years ago or will be present several million years from now. It 
can also not be excluded that a tip-RGB star was missed, for instance, because it had fallen 
outside the Spitzer coveiage, or because of an inaccuracy in the proper motion measurement, 
or for other observational reasons. To demonstrate the feebleness of any indication of a 
prematurely truncated RGB, the S-fim LF is shown again after adding a mere two stars at 
the tip of the RGB. This brings the observed LF in excellent agreement with the theoretical 
prediction. Indeed, the relative numbers of post-helium-flash red clump stars compared to 
those at similar luminosities on the RGB are well reproduced (Fig. [6]), suggesting that the 
vast majority of RGB stars do go on to undergo the helium flash. 



6. Conclusions 

We conclude that there is good agreement between the observed and expected opti- 
cal and infrared luminosity functions along the RGB in NGC6791, and that there is little 
circumstellar dust observed around RGB stars in that old metal-rich massive open cluster. 
Hence there is no direct evidence supporting the suggestion that metal-rich stars avoid the 
helium flash at the tip of the RGB and become undermassive helium-core white dwarfs as a 
result of particularly strong stellar winds. 

We thank both referees for their thorough reviews. MLB is supported by the University 
of Minnesota Louise T. Dosdall Fellowship. IM acknowledges an STFC studentship. 
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Fig. 1. — SpitzerlRAC composite of the region around NGC6791 (indicated with the cross) 
Red, yellow-green, green-cyan, and blue are 8, 5.8, 4.5, and 3.6 fim, respectively. 
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Fig. 2. — Proper mot i on m aps based on data from iMonet et al.l (l2003l ). IZacharias et al. 
(I2OO4J ) and iDias et al.l (l2006l ). which are used to select against non-members of NGC6791 
(bottom row of panels) . 
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Fig. 3. — Spitzer IRAC color-magnitude diagram of NGC6791. The RGB is well populated 
until the tip at [3.6]~ 9 mag, and the red clump is visible around [3.6]=11.5 mag. Hardly 
any star shows evidence for reddening due to a contribution of circumstellar dust emission 
at 8 fim. 
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Fig. 4. — Optical color-magnitude diagram (IStetson et al.ll2003l ). used to exclude remaining 
field stars left of the dividing line. 
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Fig. 5. — Spitzer IRAC color distribution along the RGB in NGC6791. Overplotted (dotted) 
is the same histogram, but mirrored with respect to its peak. This more clearly shows which 
stars are redder than the bulk of RGB stars. 
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Fig. 6. — Spitzer IRAC and optical V-band (jStetson et al.l l2003l ) luminosity functions in 
NGC 6791. The 8-/im luminosity function is also plotted after adding just two mo re stars near 
the tip of the RGB. A comparison with 8-Gyr Padua models (JMarigo et al.ll2008l ) shows good 
agreement, especially considering that NGC 6791 is more metal-rich. Stochastics introduced 
by sampling the small number of stars on the upper RGB can fully explain the discrepancy 
between the data and models. 



